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Sterile massive neutrinos are a natural extension of the Standard Model of elementary particles. 
The energy density of the extra sterile massive states affects cosmological measurements in an 
analogous way to that of active neutrino species. We perform here an analysis of current cosmological 
data and derive bounds on the masses of the active and the sterile neutrino states as well as on the 
number of sterile states. The so-called (3+2) models with three sub-eV active massive neutrinos 
plus two sub-eV massive sterile species is well within the 95% CL allowed regions when considering 
cosmological data only. If the two extra sterile states have thermal abundances at decoupling, Big 
Bang Nucleosynthesis bounds compromise the viability of (3+2) models. Forecasts from future 
cosmological data on the active and sterile neutrino parameters are also presented. Independent 
measurements of the neutrino mass from tritium beta decay experiments and of the Hubble constant 
could shed light on sub-eV massive sterile neutrino scenarios. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 



I. INTRODUCTION 



Solar, atmospheric, reactor, and accelerator neutrinos 
have provided compelling evidence for the existence of 
neutrino oscillations, implying non-zero neutrino masses 
(see Ref. [l| and references therein). The present data 
require the number of massive neutrinos to be equal or 
larger than two, since there are at least two mass squared 
differences (Am 2 tmos and Am 2 olar ) driving the atmo- 
spheric and solar neutrino oscillations respectively. Un- 
fortunately, oscillation experiments only provide bounds 
on the neutrino mass squared differences, i.e. they are 
not sensitive to the overall neutrino mass scale. 

Cosmology provides one of the means to tackle the ab- 
solute scale of neutrino masses. Neutrinos can leave key 
signatures in several cosmological data sets. The amount 
of primordial relativistic neutrinos changes the epoch of 
the matter- radiation equality, leaving an imprint on 
both Cosmic Microwave Background (CMB) anisotropics 
(through the so-called Integrated Sachs- Wolfe effect) and 
on structure formation, while non relativistic neutrinos in 
the recent Universe suppress the growth of matter density 
fluctuations and galaxy clustering, see Ref. 0. Cosmol- 
ogy can therefore weigh neutrinos, providing an upper 
bound on the sum of the three active neutrino masses, 
J2 m v ~ 0.58 eV at the 95% CL ||. The former bound is 
found when CMB measurements from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) are combined with 
measurements of the distribution of galaxies (SDSSII- 
BAO) and of the Hubble constant H (HST) 1 in the 
assumption of a flat universe with a cosmological con- 



1 For other recent analyses, see also Refs. 0, In- 



stant, i.e. a ACDM cosmology. 

However, the three neutrino scenario is a minimal 
scheme, and there is no fundamental symmetry in nature 
forcing a definite number of right-handed (sterile) neu- 
trino species, as those are allowed in the Standard Model 
fermion content. Indeed, cosmological probes have been 
extensively used to set bounds on the relativistic energy 
density of the universe in terms of the effective number 
of neutrinos Nf (see, for instance, Refs. [IHH ). Cur- 
rently, WMAP, SDSSII-BAO and HST data provide a 
68% CL range on Nf = 4.34^°;^ 3] in the assump- 
tion of a ACDM universe. If the effective number of 
neutrinos JV„ is larger than the Standard Model pre- 
diction of Nf = 3.046 at the Big Bang Nucleosynthesis 
(BBN) era, the relativistic degrees of freedom, and, con- 
sequently, the Hubble expansion rate will also be larger 
causing weak interactions to become uneffective earlier. 
This will lead to a larger neutron-to-proton ratio and will 
change the standard BBN predictions for light element 
abundances. Combining Deuterium and 4 He data, the 
authors of Ref. @ found Nf = 3.1±^;| at the 95% CL. 

Models with one additional ~ 1 eV massive sterile neu- 
trino, i.e. the so called (3+1) models, were introduced to 
explain LSND short baseline (SBL) antineutrino data Q 
by means of neutrino oscillations. A much better fit to 
SBL appearance data and, to a lesser extent, to dis- 
appearance data, is provided by models with two ster- 
ile neutrinos (3+2) [10|, fnjwhich can also explain both 
the MiniBooNE neutrino |12| and antineutrino data 13 1 
if CP violation is allowed (ljj . CP violation can even 
occur in (3+1) scenarios with only one relevant mass 
squared difference in the presence of non standard neu- 
trino interactions (NSI). Therefore, the (3+1) NSI model 
can also nicely explain current data [15(. While (3+1) 
and (3+2) models show some tension with BBN bounds 
on Nf, the extra sterile neutrinos do not necessarily 
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have to feature thermal abundances at decoupling. The 
first analysis of both SBL oscillation data and cosmo- 
logical data was performed by the authors of Ref. [Hj], 
where the usual full thermalization scenario for the ster- 
ile neutrino species was not assumed. Instead, the ster- 
ile abundances were computed taking into account the 
multi flavour mixing processes operating at the neutrino 
decoupling period. Robust bounds on sterile neutrino 
masses, mixings and abundances were derived. However, 
the masses of the three active neutrinos were fixe d to 
TO i "JlO- m 2 ~ V Am2 Qlar and m 3 ~ ^ Am 2 atmos . In 
Ref. [17| the authors derived the bounds on a light ster- 
ile neutrino scenario enlarging the usual thermal scenario. 
More recently, the authors of Ref. [TU have used current 
cosmological data to analyze two possible active plus ster- 
ile neutrino scenarios, one with massless active neutrinos 
(and massive steriles) and the other one with massless 
steriles states of unknown number (and massive active 
species). However, there are no cosmological bounds on 
the more natural and oscillation-data motivated scenario 
in which both the sterile and the active neutrinos have 
masses. Active neutrinos are massive; this is what oscil- 
lation data are telling us. In the same way, the LSND 
and MiniBooNE antineutrino data, if explained in terms 
of neutrino oscillations, point to the existence of massive 
sterile neutrino species. What oscillation data can not 
tell us is the absolute scale of neutrino masses and this 
is precisely what we address in this study, in the spirit of 
Ref. [T9l | . via present and future cosmological measure- 
ments. 

The paper is organized as follows. In Sec. |H]we present 
the constraints on the active and sterile neutrino masses 
and on the number of sterile species from current cosmo- 
logical data as well as from BBN measurements of light 
element abundances. Section Hill is devoted to future er- 
rors on these parameters. We describe the Fisher matrix 
method used here for forecasting errors and discuss the 
potential results from the ongoing Planck CMB mission 
combined with future BOSS and Euclid galaxy survey 
data. We also describe the induced biases on some pa- 
rameters (such as Hq and m v ) when the cosmological 
model does not account for the presence of sterile states 
to describe the data. We conclude in Sec. IIV1 



II. CURRENT CONSTRAINTS 

Here we summarize the constraints from current data 
on the active neutrino masses and on the sterile neutrino 
thermal abundance and masses. We have modified the 
Boltzmann CAMB code [2(| incorporating the extra mas- 
sive sterile neutrino parameters and extracted cosmolog- 
ical parameters from current data using a Monte Carlo 
Markov Chain (MCMC) analysis based on the publicly 
available MCMC package cosmomc[2l|. We consider here 
a flat ACDM scenario plus three {N Vs ) active (sterile) 
massive neutrino species, described by a set of cosmolog- 



ical parameters 

{L) b ,Uc,@s,T,n s ,log[10 10 A s },m v ,m Ue ,N V3 } , (1) 

where uib = flbh 2 and lu c = fl c h 2 are the physical baryon 
and cold dark matter densities, S is the ratio between 
the sound horizon and the angular diameter distance at 
decoupling, r is the optical depth, n s is the scalar spectral 
index, A s is the amplitude of the primordial spectrum 2 , 
m v is the active neutrino mass, m Ve is the sterile neu- 
trino mass and N Vs is the number of thermalized sterile 
neutrino species. We assume that both active and sterile 
neutrinos have degenerate mass spectra (m v and m Us are 
the individual masses, not the sum of the masses). The 
flat priors assumed on these cosmological parameters are 
shown in Tab. HI 



Parameter 


Prior 


Q, b h 2 


0.005-0.1 


Q. c h 2 


0.01-0.99 


Qs 


0.5-10 


T 


0.01-0.8 


n s 


0.5-1.5 


ln(10 10 A s ) 


2.7-4 


m„ s [eV] 


0-3 


m v [eV] 


0-3 




0-6 



TABLE I: Flat priors for the cosmological parameters consid- 
ered here. 



Our basic data set is the seven-year WMAP data 
3, 22] (temperature and polarization) with the routine 
for computing the likelihood supplied by the WMAP 
team. We consider two cases: we first analyze the 
WMAP data together with the luminous red galaxy 
clustering results from SDSSII (Sloan Digital Sky Sur- 
vey) [23d and with a prior on the Hubble constant from 
HST (Hubble Space Telescope) |24| . referring to it as 
the "runl" case. We then include to these data sets Su- 
pernova la Union Compilation 2 data [25j, and we will 
refer to this case as "run2". In addition, we also add 
to the previous data sets the BBN measurements of the 
4 He abundance, considering separately helium fractions 
of Yp 1 = 0.2561 ± 0.0108 (see Ref. (2ft) and of Y3 = 
0.2565 ± 0.0010 (stat.) ±0.0050 (syst.) from Ref. [27| . 
Finally, we also consider the Deuterium abundance mea- 
surements log(D/H) = 4.56 ± 0.04 from Ref. H<|. 

It is important to clarify that CMB anisotropics also 
depend on the value of Y v but since Y p is constrained 
loosely by current CMB/LSS data, it is consistent to fix 
it to value K, = 0.24 in the CMB runs and to consider 



2 The pivot scale assumed in this study corresponds to kg = 
0.05 Mpc" 1 . 
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it as an independent parameter constrained by BBN ob- 
servations. 

Given a cosmological model, we predict the theoretical 
primordial abundance of Y p and \og(D/H) by making 
use of the public available PArthENoPE BBN code (see 
30]). 

Since running cosmomc and getting at the same time 
the theoretical predictions from Parthenope for BBN 
would be be exceedingly time-consuming we perform im- 
portance sampling obtaining the predicted values for Y p 
and \og(D/H) with an interpolation routine using a grid 
of Parthenope predictions for each (u>b, N Vs ), as in |28| . 

Parameter 68% CL(rl) 95% CL(rl) 68% CL (r2) 95% CL (r2) 



m v [eV] 
m Vs [eV] 



< 2.5 

< 0.13 

< 0.22 



< 4.1 

< 0.30 

< 0.46 



< 2.0 

< 0.10 

< 0.20 



< 3.2 

< 0.20 

< 0.50 



TABLE II: ID marginalized bounds on the active and sterile 
neutrino parameters using the two combinations of data sets 
described in the text (rl refers to "run 1" and r2 refers to 
"run 2", respectively). 



[26] y p 2 [27] Kp 1 + D [29] Y% + D [29] 

N Va < 2.3 < 1.7 < 1.7 < 1.4 

m v [eV] < 0.17 < 0.15 < 0.15 < 0.15 

m„ 3 [eV] < 0.62 < 0.67 < 0.69 < 0.68 

TABLE III: ID marginalized 95% CL bounds on N u , , m v , 
and m„ after combining the results of "run 2" with those 
coming from different measurements of BBN light element 
abundances. 



Table HH shows the ID marginalized bounds on N v , 
m Vs and m v arising from the two different analyses per- 
formed here on cosmological data sets. The marginalized 
limits have been computed setting a lower limit of in all 
the three neutrino parameters here explored. The bounds 
obtained on the parameters associated to the dark mat- 
ter candidates considered here are consistent with those 



obtained in Ref. [31| after taking into account the differ- 
ences in the thermal abundances of sterile neutrinos and 
QCD thermal axions. When we marginalize over all the 
cosmological parameters, see Tab. HH the 95% CL upper 
bound for N v , is 4.1 (3.2) using "runl" ("run2") data 
sets. Therefore, current cosmological data does not ex- 
clude at the 95% CL the existence of ~ 2 sterile neutrino 
species with sub-eV masses plus three sub-eV active mas- 
sive neutrinos. It would be interesting to further explore 
if a model with sterile neutrinos is preferred over the 
model with only three active neutrinos. The results here 
are also in very good agreement with those of Ref. [l8[ 
even if in the former analysis the two species, i.e. the 
active and sterile neutrino states, were not considered to 
be massive at the same time. 



Table Mil shows the 95% ID marginalized bounds on 
N Vs , m„ s and m„ arising when different combinations of 
BBN light clement abundances measurements are com- 
bined with "run 2" results. Note that when measure- 
ments of the 4 He abundance are added to CMB, galaxy 
clustering and SNIa data, the 95% CL upper limit on 
N v , is 2.3 (1.7) if Y p = 0.2561 ± 0.0108 (Y p 2 = 0.2565 ± 
0.0010 ± 0.0050) is assumed. Since the number of sterile 
species after adding BBN constraints is smaller than be- 
fore, the sterile (active) neutrino masses can get slightly 
larger (smaller) values, since BBN data is insensitive to 
the dark matter density in the form of massive neutri- 
nos at late times. The combination of Helium and Deu- 
terium abundance measurements compromises the via- 
bility of (3+2) models, leading to N Vs < 1.7 - 1.4 at 
the 95% CL. However, the two sterile states might not 
have thermal properties at decoup ling and evade BBN 
constraints. A complete analysis [32j including sterile 
neutrino mixing parameters and recent reactor neutrino 
oscillation results [33| is mandatory. 

Figure [TJ top panel, depicts the 68% and 95% CL al- 
lowed contours in the m u -N v ^ plane. The blue (red) 
contours denote the allowed regions by "runl" ("run2") 
data sets. Notice that there exists a degeneracy between 
these two quantities. This degeneracy is similar to the 
one found by the authors of Ref. [18| . When the mass en- 
ergy density in the form of massive neutrinos is increased, 
the number of extra relativistic species must also be in- 
creased to compensate the effect. This will be the case for 
massless sterile species. In our analysis, the degeneracy 
is milder since sterile neutrinos are massive and therefore 
they behave as an additional dark matter component at 
late times. The degeneracy will show up when the active 
neutrinos have relatively large masses, since, in that case, 
a tiny amount of sterile neutrino masses will be allowed. 
The sterile states will then behave as relativistic particles 
at the decoupling era and will compensate the effect of a 
large active neutrino mass. 

Figure [U middle panel, depicts the 68% and 95% CL 
allowed contours in the m„-m„ s plane. There exists a 
very strong anticorrelation between these two quantities, 
since both contribute to the dark matter energy density 
at late times and therefore if the mass of the sterile neu- 
trino states grows, the mass of the active ones must de- 
crease. The situation is analogous to that of QCD ther- 
mal axions and massive (active) neutrinos, see Ref. [31j ]. 

The bottom panel of Fig. Q] depicts the 68% and 95% 
CL allowed contours in the N Vs -m Vs plane. In this case, 
the larger the sterile neutrino mass is, the lower its ther- 
mal abundance must be, as expected. 



III. FUTURE CONSTRAINTS 

We present here the constraints on the neutrino sector 
parameters explored in this work from future CMB and 
galaxy survey measurements, making use of the Fisher 
matrix formalism, see also Ref. 34] for a recent analy- 
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FIG. 1: The top, middle and bottom panels show the 68% and 
95% CL constraints on the plane m v -N V3 , m v -m v ^ and m Vs - 
N„ s , respectively. The blue (red) contours denote the allowed 
regions by "runl" ( "run2" ) data sets, see text for details. The 
masses of the sterile and active neutrinos are both in eV units. 



sis. We also compute the potential shifts in the different 
cosmological parameters when the sterile neutrino pa- 
rameters are neglected in the analysis. 



A. Methodology 

The Fisher matrix is defined as the expectation value 
of the second derivative of the likelihood surface about 
the maximum. As long as the posterior distribution for 
the parameters is well approximated by a multivariate 
Gaussian function, its elements are given by (35rl37| 



F a p = ^Tr [C-^cC-iCj] , 



(2) 



where C = S + N is the total covariance which consists 
of signal S and noise N terms. The commas in Eq. (j5J) 
denote derivatives with respect to the cosmological pa- 
rameters within the assumed fiducial cosmology. Our 
fiducial model is a ACDM cosmology with five parame- 
ters: the physical baryon and CDM densities, a;& = Slbh 2 
and uj c — tt c h 2 , the scalar spectral index, n s , h (being 
the Hubble constant Hq = 100 h km Mpc _1 s~ 1 ) and 
the dimensionless amplitude of the primordial curvature 
perturbations, A s (see Tab. IIVI for their values). Fur- 
thermore, we add to the ACDM fiducial cosmology three 
additional parameters for the neutrino sector: the mass 
of active neutrinos m„ , the mass of sterile neutrinos m Vs 
and the number of sterile neutrino species N Vs . Notice 
that, for simplicity, we have kept fixed the reionization 
optical depth r since it has no impact on large scale struc- 
ture data and we do not expect a strong degeneracy be- 
tween t and the neutrino parameters, see Ref. [38j |. We 
assume that both active and sterile neutrinos have a de- 
generate spectrum and that the sterile species are fully 
thermalized. The fiducial values of the neutrino param- 
eters are listed as well in Tab. IIVI and they are based on 
the constraints from current data presented in the previ- 
ous section, from which we conclude that m„ = 0.1 eV, 
m Ut < 0.5 and N Us = 1, 2 are within the allowed regions 
for these parameters. 



tl b h 2 


n c h 2 


n s 


h 


A s 


m v [eV] 


m„ s [eV] 




0.02267 


0.1131 


0.96 


0.705 


2.64 • 10" 9 


0.1 


0.1-0.5 


1-2 



TABLE IV: Values of the parameters in the fiducial models 
explored in this study. 

We compute the CMB Fisher matrix to obtain fore- 
casts for the Planck satellite (3{|. We follow here the 
method of Ref. |io| . considering the likelihood function 
for a realistic experiment with partial sky coverage, and 
noisy data 



-21n£ = ^(2£+l) /.. 



rBB , / 

I sky m 



f-iBB 



+ 



(Cj E f 



+ 



/ fTT cEE 
y ' sky J sky 



TT f~\EE 



~<TTfiEE 



QTTqEE 



( C I E ) 2 



+ 



. rBB W 
TJ sky g 



fTT fEE 
J sky J sky 



rBB 
J sky 



(3) 



and compute its second derivatives to obtain the corre- 
sponding Fisher matrix 



rCMB 



d 2 L 

dpadpp 



(4) 



In Eq. © Cf Y = Cf Y + Af e XY with Cf Y the tempera- 
ture and polarization power spectra (X, Y = {T, E, B}) 
and Mi the noise bias. Finally, f~jZ is the fraction of 
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observed sky which can be different for the T-, E-, and 
-B-modes. 

For the galaxy redshift survey Fisher matrix, we follow 
the prescription of Ref. [4l|. Assuming the likelihood 
function for the band powers of a galaxy redshift survey 
to be Gaussian, the Fisher matrix can be approximated 



? LSS 



ainP gg (fc)91nP gg (fc) 



* d\nP m 



dpf3 



V cS (k) 



dk 



2(2tt) s 



(5) 



-1 J k„ 
2, 



8p a 



dp/3 



V eS (k,iJ,) 



2nk 2 dkdfi 



2(2tt) 3 ' 

where V e f{ is the effective volume of the survey: 

nP(k, n) 



V eS (k,fj.) 



1 2 



Vs, 



(6) 



_nP(k,fi) + 1 

/i being the cosine of the angle between the vector along 
the line of sight and k and n being the galaxy number 
density, which is assumed to be constant throughout the 
survey. The linear redshift-space galaxy power spectrum 
P gg is related to the real-space linear power dark matter 
spectrum Pdm as 

Pg g (fc)-Pdm(fc)(fe+//i 2 ) 2 (7) 

where b is the bias relating galaxy to dark matter over- 
densities in real space and / is the linear growth factor. 
Both the bias and the growth factor are assumed to vary 
in each redshift bin and are considered as additional pa- 
rameters in the Fisher analysis of galaxy survey data. 

We consider here two redshift surveys: the BOSS 
(Baryon Oscillation Spectroscopic Survey) [42j and the 
Euclid [43|, |3] experiments. For the BOSS survey we 
assume a sky area of 10000 deg 2 , a redshift range of 
0.15 < z < 0.65 and a mean galaxy density of 2.66 x 10~ 4 . 
For Euclid we consider an area of 20000 deg 2 , a redshift 
range of 0.15 < z < 1.95 and a mean galaxy density 
of 1.56 x 10~ 3 . We divide the surveys in redshift bins 
of width Az = 0.1 (a value that is much larger than 
standard redshift spectroscopic errors), set /c max to be 
0.1/i/Mpc and fc m j n to be greater than 2n/ AV 1 ^ 3 , where 
AV is the volume of the redshift shell. 

Combining the Planck and redshift survey Fisher ma- 
trices (F a p = F^ s + F^ B ) we get the joint constraints 
for f2;,/i 2 , fl c h 2 , n s , Hq, A s , m„, m Us and N Us , after 
marginalizing over the bias b and the growth factor /. 
The 1-ct error on parameter p a marginalized over the 
other parameters is a(p a ) = \J {F~ v ) aa , F^ 1 being the 
inverse of the Fisher matrix. 



B. Results 

Tables El and EH contain the l -a marginalized fore- 
casted errors on the cosmological parameters for a fidu- 
cial cosmology with m v = 0.1 eV, m Us — 0.3 eV and 



N Vs — 1 and 2, respectively. We illustrate the results 
of our Fisher analysis for both BOSS and Euclid galaxy 
redshift survey data combined with Planck CMB mea- 
surements. Note that the errors on the pure ACDM 
model parameters, i.e ttbh 2 7 Q c h 2 , n s , h and A s are al- 
ways around or below the 1% level. The error on the 
active neutrino mass is around 60% for BOSS and half 
for Euclid. The error on the number of sterile neutrino 
species is always smaller than 25%. Regarding the error 
on the sterile neutrino mass, it can reach 100% relative 
errors for BOSS plus Planck data. Naively, one would 
expect that the BOSS and Euclid errors are related by 
a factor of ^/Vboss /^Euclid (being V the volume of the 
survey) when the shot noise is subdominant. However, 
in practice, the forecasted errors on the pure ACDM pa- 
rameters are sometimes similar for the BOSS and Euclid 
cases, which implies that those parameters are mainly de- 
termined by CMB measurements. Of course this is not 
the case for the active and sterile neutrino masses, whose 
errors are mainly driven by galaxy clustering data and 
differ by a factor of ~ 2 — 3 for BOSS and Euclid cases, 
as naively expected. A word of caution is needed here: 
while computing the errors on the active and sterile neu- 
trino masses and on the sterile neutrino abundances, a 
ACDM scenario has been chosen as fiducial cosmology. 
These errors can change if the equation of state of the 
dark energy component is allowed to vary (45j and/or 
interactions between the dark matter and dark energy 
sectors are stwiched on WO, U3I . 



We also present here the joint constraints in a two- 
parameter subspace (marginalized over all other cos- 
mological parameters) to study the covariance between 
the sterile neutrino masses and/or abundances and the 
other cosmological parameters considered in this work. 
We have explored several possible scenarios with differ- 
ent sterile neutrino masses and thermal abundances (see 
Tab. IIV[) . However, for the sake of simplicity, we illus- 
trate here only the case N l/s = 1, m Vs = 0.3 eV and 
m v = 0.1 eV. 

Figure [21 left panel, shows the correlation between the 
number of sterile species N Va and the active neutrino 
mass m„. The expected error on the number of sterile 
species is very similar for BOSS and Euclid data, which 
indicates that the constraints on N Ug arise mostly from 
Planck CMB measurements. Since the total energy den- 
sity in the form of massive neutrinos is the sum of the 
active plus sterile contributions, a higher neutrino mass 
is compensated with a lower abundance of massive sterile 
species. The 1— a marginalized error on N Vg from Planck 
plus BOSS (Euclid) data is 0.26 (0.1), see 3 Tab. El The 
right panel of Fig. [2] shows the correlation between the 
masses of sterile and active neutrinos. As expected from 
the results presented in Fig. Q] (middle panel) and as pre- 
viously explained, higher active neutrino masses are al- 
lowed for very low values of the sterile neutrino masses. 
The 1-a marginalized errors on the massive species m Vs 
and m„ from Planck plus BOSS (Euclid) data are 0.25 
(0.08) eV and 0.06 (0.03) eV respectively. If nature has 
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chosen an active neutrino with mass ~ 0.1 eV, BOSS (Eu- 
clid) data, combined with CMB Planck measurements, 
could provide a 1.5-c (3-ct) detection, even in the pres- 
ence of massive sterile species. 



Parameter BOSS+PLANCK EUCLID+PLANCK 



n b h 2 


0.7% 


0.3% 


tt c h 2 


2.9% 


1.3% 


ln(10 10 A s ) 


0.7% 


0.4% 


h [km/s/Mpc] 


1.4% 


0.7% 


n s 


0.6% 


0.3% 


m v [eV] 


63.1% 


28.0% 


m Vs [eV] 


83.2% 


26.2% 




25.9% 


10.6% 




m y (eV) m v (eV) 



FIG. 2: The empty (filled) contours denote the 68%, 95% 
and 99.73% CL regions for Planck plus BOSS (Euclid) data. 
The neutrino parameters in the fiducial model are jV„ 3 = 1, 
m Va = 0.3 eV and m v —0.1 eV. 



TABLE V: lcr marginalized relative errors for a fiducial cos- 
mology with Nv s =l, m v =0.1 eV and m Vs —0.3 eV. 



Parameter BOSS+PLANCK EUCLID+PLANCK 



n b h 2 


0.7% 


0.3% 


tt c h 2 


1.5% 


1.7% 


ln(10 10 A) 


0.4% 


0.4% 


h [km/s/Mpc] 


1.4% 


0.8% 


n s 


0.5% 


0.4% 


m v [eV] 


64.9% 


35.9% 


m„ 3 [eV] 


41.9% 


16.4% 


N Va 


10.2% 


7.5% 



TABLE VI: l—o marginalized relative errors for all parameters 
for a N Vs =2, m„=0.1 eV and m„ 3 =0.3 eV fiducial cosmology. 




m, (eV) N„ 



FIG. 3: The empty (filled) contours denote the 68%, 95% and 
and 99.73% CL regions for Planck plus BOSS (Euclid) data. 
The neutrino parameters in the fiducial model are N Vs — 1, 
m„ 3 = 0.3 eV and m„ = 0.1 eV. 



C. Cosmological parameter shifts 



Figure [31 left panel, shows the correlation between the 
active neutrino mass m„ and the cold dark matter en- 
ergy density il c h 2 . Notice that the extraction of the 
cold dark matter component arise mostly from Planck 
CMB measurements. At late times, neutrinos contribute 
as an additional ingredient to the dark matter fluid and 
therefore a higher neutrino mass is compensated by a 
lower cold dark matter energy density. The right panel 
of Fig. [3] shows the correlation between cold dark mat- 
ter and the sterile neutrino abundance. These two pa- 
rameters are mostly extracted from CMB Planck data 3 . 
The sterile neutrinos considered here with 0.3 eV masses 
are relativistic at decoupling. A higher number of rel- 
ativistic species will shift to a later period the matter 
radiation equality era and also enhance the first CMB 
acoustic peak. These effects can be compensated with a 
higher cold dark matter energy density, as shown by the 
positive correlation among the two parameters. 



However, the addition of galaxy clustering measurements help in 
breaking degeneracies 



In order to test the capabilities of future experiments 
to discriminate between different theoretical models, re- 
gardless of their parameters, we follow here the method 
of Ref. 0. 

The idea is the following: if the data is fitted assuming 
a model M ' with n' parameters, but the true underlying 
cosmology is a model M characterized by n parameters 
(with n > n' and the parameter space of M including 
the model M' as a subset), the inferred values of the n' 
parameters will be shifted from their true values to com- 
pensate for the fact that the model used to fit the data is 
wrong. In the case illustrated here, M will be the model 
with massive sterile neutrinos and M' the one without 
massive sterile neutrinos. While the first n' parameters 
are the same for both models, the remaining n — n' = p 
parameters in the enlarged model M are accounting for 
the presence of massive sterile neutrinos, i.e. m Vs and 
N Vs . Assuming a gaussian likelihood, the shifts of the 
remaining n! parameters is given by [48(: 

59' a = -(F'-%pG K 5ip c a, = 1 . . .n', ( = n'+l ...n, 

where F' represents the Fisher sub-matrix for the model 
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M' without massive steriles and G denotes the Fisher 
matrix for the model M with m Vs , N Ve > 0. 

We have computed the shifts induced in the cosmologi- 
cal parameters in several true cosmologies with a number 
of sterile neutrinos N Ua — 1,2 of masses m„ s = 0.1,0.3 
and 0.5 eV. The mass of the active neutrino has been 
kept to 0.1 eV. These cosmologies are then wrongly fitted 
to a cosmology without sterile massive neutrino species. 
While certain parameters are exclusively measured by 
CMB probes or by the combination of CMB and other 
cosmological data sets (like £l c h 2 , £lbh 2 , n s and A s ), there 
are other parameters such as the Hubble constant ffo or 
the active neutrino mass m v which can be determined by 
other experiments. Then it is possible to verify the cos- 
mological model assumptions by comparing the values of 
Hq and m v extracted from CMB and LSS cosmological 
data to the values of these parameters obtained by other 
experiments, as missions devoted to measure the Hub- 
ble constant and tritium beta decay experiments 4 . The 
former experiments measure the electron neutrino mass 
m„ e , which, in practice, when considering three active 
massive neutrinos, reads: 

< = E I^K . (9) 

i=l,3 

being U e i the first-row elements of the Pontecorvo-Maki- 
Nakagawa-Sakata leptonic mixing matrix. In the case 
of additional 7V„ s massive sterile neutrino species, m„ e 
would be given by 

< = E \ u *\ m * , (io) 

Given the current best-fit values of the sterile-electron 
neutrino mixing terms \U es \ < 10 _1 [l6| and the sub-eV 
sterile neutrino masses considered here, we neglect the 
contribution of the sterile neutrino species to m Ve . In 
the following, we apply the usual constraints on m„ c in 
our cosmological scenarios even if they contain massive 
sterile neutrino species. We therefore neglect the capa- 
bility of beta decay experiments of measuring the indi- 
vidual neutrino masses and mixings. For a recent study 
of the KATRIN potential for sterile neutrino detection, 
see Ref. 

For sterile neutrino masses m„ s ~ 0.5 eV and N Vs = 
1,2, the shifts induced in Hq are very large, for both 
BOSS and Euclid experiments combined with CMB 
Planck data. The reconstructed value of Hq is within the 
range ^ 20 — 50 km/s/Mpc, values which are in strong 
disagreement with current measurements of the Hubble 



4 Neutrinoless double beta decay provides also a bound on the so- 
called effective neutrino mass (m) = | U^rrii]. However, these 
bounds apply only in the case that neutrinos have a Majorana 
nature. Therefore, we focus on tritium beta decay constraints 
which apply regardless of the Dirac vs Majorana nature of the 
neutrino. 



parameter from HST [24], 50] . The reconstructed value of 
the active neutrino mass is also in some cases m v ~ 2 eV 
which is the current 95% CL limit from tritium /3-decay 
experiments [5lj . Consequently, after combining near fu- 
ture BOSS and Planck data one would conclude that 
the cosmological model assumed with m Vs ~ 0.5 eV and 
N Va = 1,2 is wrong. The same situation will arise when 
m„ s ~ 0.3 eV and two sterile massive species, N„ s = 2. 



Parameter 


Fiducial Reconstructed Shift (%) 


H [km/s/Mpc] 
m„ [eV] 


70.5 50.5 28% 
0.30 0.98 230% 



TABLE VII: Shifted values and relative changes for the pa- 
rameters Ho and ra„ when the true cosmology has N Us = 1, 
m Vs = 0.3 eV and m„ = 0.1 eV but BOSS plus Planck data 
are fitted to a cosmology with no sterile massive neutrino 
species. 



Parameter 


Fiducial Reconstructed Shift (%) 


H [km/s/Mpc] 
m v [eV] 


70.5 65.0 8% 
0.30 0.48 60% 



TABLE VIII: Shifted values and relative changes for the pa- 
rameters Ho and m v when the true cosmology has N Us — 1, 
m„ 3 = 0.3 eV and m„ =0.1 eV but Euclid plus Planck data 
are fitted to a cosmology with no sterile massive neutrino 
species. 

For m„ 3 ~ 0.3 eV and N Vs = 1, the shifts using both 
BOSS and Euclid data are reported in Tabs. IVIII and 
IVIIII While the shift induced in the Hubble constant is 
very large for the BOSS case, for Euclid that shift is 
still consistent with current estimates of Hq. A number 
of experiments (HST, Spitzer, GAIA and JWST HI) 
are expected to measure -ffo with 2% uncertainty in the 
next decade and an inconsistency between the inferred 
ffo values from these experiments and those from the 
cosmological probes considered here could point to the 
existence of additional sterile neutrino species. On the 
other hand, the aim of the tritium beta decay experi- 
ment KATRIN 53] is a sensitivity of m Ve . < 0.2 eV at 
90% CL in case of a null result or a 5ct discovery potential 
for m„ e > 0.35 eV. Therefore, the reconstructed values of 
m v = 0.48 eV (Euclid plus Planck) and 0.98 eV (BOSS 
plus Planck) could be easily testable by the KATRIN ex- 
periment. Similar results are obtained for smaller sterile 
neutrino masses m Ug ~ 0.1 eV with a higher number of 
sterile species N Us = 2. 

For smaller sterile neutrino masses m Vs ~ 0.1 eV and 
N Us = 1, the shift induced in ffo is larger than 2% for 
both BOSS and Euclid data (combined with Planck). 
Therefore it would still be possible to check the fiducial 
cosmology with future measurements of Hq. The shift 
induced on the active neutrino mass using Euclid data 
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is negligible and this means that it would be possible to 
recover the true value of the active neutrino mass even 
if the data is fitted to the wrong cosmology. Thus, the 
combination of Planck and Euclid data would not lead to 
an inconsistency between active neutrino mass estimates 
from Planck and Euclid on the one hand, and beta de- 
cay experiments on the other hand. Regarding BOSS 
plus Planck data however, the shift induced in the ac- 
tive neutrino mass is of the order of 100% and the 
comparison with an independent measurement of m v as 
that performed by KATRIN could test the validity of the 
cosmological model assumptions. 

We have shown above that if the true N Ve = 1,2, 
wrongfully assuming N Vs — would lead to discrep- 
ancies between the cosmological probes considered here 
(large scale structure and CMB) and independent mea- 
surements of Hq and m v . Of course, another clear indi- 
cator that the assumed model is incorrect is simply that 
the N Vs = would likely provide a bad fit to the large 
scale structure and CMB data themselves. However, the 
induced bias discussed above would provide a useful ex- 
tra check when independent measurements of Hq and/or 
m v are available. In addition, the bias calculation shows 
that even if one is not interested in the sterile neutrinos 
per se, not taking them into account could lead to very 
wrong conclusions about the other cosmological parame- 
ters. 



IV. SUMMARY 

Neutrino oscillation experiments have brought to light 
the first departure from the Standard Model of particle 
physics, indicating that neutrinos have non zero masses 
and opening the possibility for a number of extra ster- 
ile neutrinos. LSND and MiniBooNE antineutrino data 
require these extra sterile species to be massive. Much 
effort has been devoted in the literature to constrain the 



so called (3+1) (three active plus one sterile) and (3+2) 
(three active plus two sterile) models. 

Cosmology can set bounds on both the active and ster- 
ile neutrino masses as well as on the number of ster- 
ile neutrino species. We have explored here the current 
constraints on these parameters in the most natural sce- 
nario which corresponds to the case in which both the 
active and sterile neutrinos are massive particles. We 
find that models with two massive sub-eV sterile neu- 
trinos plus three sub-eV active states are perfectly al- 
lowed at the 95% CL by current Cosmic Microwave Back- 
ground, galaxy clustering and Supernovae la data. The 
bounds derived here were obtained in the context of a 
ACDM cosmology and other scenarios with a dark en- 
ergy component could allow for larger neutrino masses 
and/or abundances. We have also shown that Big Bang 
Nucleosynthesis Helium-4 and deuterium abundances ex- 
clude (3+2) models at the 95% CL. However, the extra 
sterile states do not necessarily need to feature thermal 
abundances at decoupling. Their precise abundances are 
related to their mixings with the active neutrinos in the 
early universe. 

We have also forecasted the errors on the active and 
sterile neutrino parameters from Planck and galaxy sur- 
vey data. Future cosmological data are expected to 
measure sub-eV active and sterile neutrino masses and 
sterile abundances with 10 — 30% precision, for sub-eV 
(0.5 eV> m Vs > 0.1 eV) sterile neutrino masses. We 
have also shown that the presence of massive sterile neu- 
trinos in the universe could be inferred from inconsisten- 
cies among the values of Hq obtained from cosmic mi- 
crowave and galaxy clustering probes and those arising 
from independent measurements of the Hubble constant 
over the next decade. The validity of the cosmological 
assumptions could also be tested by comparing cosmo- 
logical measurements of the active neutrino mass with 
those obtained from tritium beta decay experiments. 
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